The multilevel pairing Hamiltonian versus the degenerate case 
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Abstract 

We study the pairing Hamiltonian in a set of non degenerate levels. First, we review in the 
path integral framework the spontaneous breaking of the U(l) symmetry occurring in such a 



system for the degenerate situation. Then the behaviors with the coupling constant of the ground 
state energy in the multilevel and in the degenerate case are compared. Next we discuss, in the 
multilevel case, an exact strong coupling expansion for the ground state energy which introduces 
the moments of the single particle level distribution. The domain of validity of the expansion, 
which is known in the macroscopic limit, is explored for finite systems and its implications for 
the energy of the latter is discussed. Finally the seniority and Gaudin excitations of the pairing 
Hamiltonian are addressed and shown to display the same gap in leading order. 
PACS: 21.60.-n, 21.30 Fe, 24.10.Cn; Keywords: Pairing interaction; Bosonization. 

1. Introduction 

The pairing Hamiltonian 
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31 

j m=-j i,j 



# = E E ejX^Xjrn- g^xm^Uh (1) 



where (note that this definition differs from the one of Refs. [l|2|3j by a normalization factor 

l/VO") 



Ai — "Tfr ^ m Xj- m \j, m , (2) 

V lL i m>0 

represents a simplified version of the BCS model of superconductivity since it assumes a constant 
two-body matrix element of the pairing force between all the single particle levels. 
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While BCS is generally applied to a macroscopic system, the Hamiltonian (pQ) is also suitable 
for dealing with finite systems. Indeed it has been successfully applied to atomic nuclei where N 
fermions live in L single particle levels, each with pair degeneracy JL- = j + 1/2. In this case in 
dH) AT, A are the usual creation and annihilation fermion operators, m is the third component of 
the angular momentum j, the ej are the single-particle energies (assumed to be m- independent) 
and g is the strength of the pairing force. 

One recovers BCS from ([1]) by taking Qj = 1 for all the L levels and letting L and N become 
very large (as appropriate to a macroscopic system, e.g. electrons in a metal), but keeping their 
ratio constant and gL finite. 

As is well-known in this limit the ground state of the grandcanonical version of ([T]) reads H] 

\BCS>= [] ( Ufc + Vfc At A f _ a )|0> , (3) 

fc>0 

|0 > being the particle vacuum and Uk, Vk parameters taken to be real. The operators in (J3j) 
create a pair of particles in time reversal states. 

The state (|3|) is not invariant under the group U(l) of the global gauge symmetry respected 
by (PQ), hence it is associated with the spontaneous breaking of the particle number conservation. 

When applied to a system of a finite, fixed number ./V of particles ([3]) becomes an approxi- 
mation, valid, however, to the order 1/N. This follows from the finding of Cambiaggio et al. [5] 
that the Hamiltonian (fTJ) corresponds to an integrable system: hence it is solvable (classically 
by quadrature, according to the Liouville theorem). Indeed Richardson [6] was able to find out 
the system of algebraic equations yielding the exact solution which, for a system made up of 
n = N/2 (N even) pairs of fermions living in L single particle levels, reads [2j 



\n >= C n J] Bt\0 > (4) 
being 

= E 2e^W Aj ■ (5) 

j = l ZK 7 

In ([!]) C n is a normalization constant and in ([5]) the quantities E v (y = 1 • • • n) are solutions of 
the following system of n equations (the Richardson equations) 

1 L O n 2 

- = Y - 1 V - (6) 

whereas the Aj (defined in ) are referred to as hard-core boson operators and satisfy the 
commutation relations 



Aj,Aj, 



V (l-$ , (7) 
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being Nj = J2m=-j ^jm^jm the particle number operator of the j-level. 

In Richardson's framework, the eigenvalues £ of (H|) (in particular the ground state energy) 
are obtained according to 

n 

£(n) = £>„. (8) 

v=l 

A question arises on the relationship between and ([3]) in the limit of N and L very large, 
the only situation where, strictly speaking, the phase transition associated with the spontaneous 
symmetry breaking occurs. We shall later comment on this issue, already addressed by Roman et 
al. [7] who, in analyzing this limit, showed that the system of equations ^ becomes a non linear 
integral equation and identified the order parameter associated with the spontaneous breaking 
of particle number conservation. This turns out to be the BCS gap, namely 

A = 9 U P V » ■ ( 9 ) 

Ai>0 

In general ^ can only be dealt with numerically. However it is trivial to solve it when only 
one single particle level of energy e and degeneracy ft is available to the fermions. In fact the 
system then reduces to the form 

1 ft A 2 

9 2e ~ E - M( ^ )=1 E » ~ E » 
or 

E u = 2e-gn-2g £ + Age £ , (11) 

which, upon summation over the index v, immediately yields the well-known expression for the 
ground state energy 

n 

8 = ^2 E v = 2en - gn (O - n + 1) . (12) 

Note that in the above the term quadratic in n has a positive sign, reflecting the action of the 
Pauli principle. 

In the L = 1 case it is also easy to get the excitation spectrum of the Hamiltonian ([1]) which 
can only correspond to the breaking of pairs. Indeed this occurrence can be formally accounted 
for by changing ft into ft — 2s and n into n — s, s (the so-called pair seniority) counting the 
number of broken pairs. One ends up with the remarkably symmetric expression 

n 

£ s = £ E v = 2en - gn (ft - n + 1) + gs (ft - s + 1) , (13) 

which yields 2en, as it should, for s = n, namely when all pairs are broken, thus becoming blind 
to the pairing interaction. 
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2. The path integral approach in the degenerate case 

The equations ©, while exact, do not transparently convey the occurrence of the spontaneous 
symmetry breaking taking place in the system. 

To shed light on this point it turns out that a path integral approach, while difficult, is more 
suited. This scheme has been followed in Ref. [1] starting from the discretized Euclidean action 
corresponding to the Hamiltonian ([1]), which in the L = 1 case reads [8] 

S = t J2 \-gnA(t)A(t-l)+ [Am(i)(v t + + e)A m (i-l)] \, (14) 

t=-N /2 { tn=-j J 

r being the time spacing, Nq the number of points on the time lattice, the energy e is mea- 
sured with respect to the chemical potential to select a sector of N fermions and A, A are odd 
Grassmann variables. In terms of the latter 

A = ~7^J2 (-l) j-m A-mA m , (15) 
v^ m>0 

V"m>0 

and the discretized time derivative is defined as follows 

(v±/)(t) = ±~[/(i±l)-/(t)] . (17) 

Next the action (j!4[) is dealt with through the Hubbard-Stratonovitch transformation by intro- 
ducing a bosonic field r\ cast in the polar representation: 

r, = ^Tpe 2W , fj = ^e~ 2W . (18) 



In the end the following effective action, equivalent to (I14p and describing the same physics of 
the Hamiltonian ([1]), 

S cS = r^ 5 p-Tr In (-q~ q + + g 2 p) (19) 

t 

is obtained. In (|19[) the trace is taken both on the time and on the quantum number m and 
g ± = e* ie Vfe ±ie ± e . (20) 
Now in the above we see that the field 9 

a) appears in the action only under derivative, 

b) lives in the coset space of the broken group 17(1) with respect to the unbroken group Z<z- 
This follows from the requirement of making the change of variable (|18|) to be one to one. 
Indeed for this to occur it must be < 6 < it. 
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c) Moreover, at variance with the Hamiltonian (pQ), the U(l) symmetry is now non-linear ly 
realized in the invariance of the action (|19p under the substitution 6 —* 6 + a , a being 
time independent. 

The above features characterize a Goldstone field (see [§]). 

These findings show that in the degenerate case the system lives in a regime where the U(l) 
symmetry is spontaneously broken not only at the macroscopic, but at the microscopic scale as 
well. 

According to the Goldstone theorem when a continuous symmetry is broken a bosonic field 
appears with a dispersion relation which vanishes in the thermodynamic limit. 

Of course in a finite system, like a deformed nucleus, where the 0(3) symmetry is sponta- 
neously broken , the Goldstone boson shows up as a rotational band whose frequency is not 
vanishing, but, however, tends to zero as the system becomes larger and larger. Likewise in 
the pairing problem the £7(1) symmetry is broken and in the degenerate case the energy of the 
associated Goldstone mode depends upon the pair number n according to (fT2|) . Actually setting 
(see P) 

v = n- [n ] , (21) 

where the square bracket means integer part and 

^+1 e , , 

no = -= 22 

2 g 

is the value of n where f)12|) reaches its minimum, one sees that the energy of the Goldstone 
mode goes like g, whereas the energy associated with the variable s, hence corresponding to 
the breaking of pairs, goes like g£l, thus lying at a much higher energy. Microscopically the 
Goldstone mode corresponds to the addition or to the removal of a pair of fermions. 

In the path integral framework another point is of relevance. In dealing with S e s with the 
saddle point expansion one starts by searching for its minimum at constant fields. It turns then 
out that the minimum of the effective action occurs for p = p, being [10] 



P 



I i A 2 

i r, ,_ 2 i ^ 



(gny-4e 2 = — , (23) 

g 2 



(29) 2 

where A, see ©, is the well-known gap characterizing the BCS theory in the L = 1 case. Thus 
p, but for a factor which renders it dimensionless, coincides with the gap. Furthermore, choosing 
e = 0, <|23D yields gQ = 2A. 

3. The Richardson approach in the multilevel case 



We now study the case when L single particle levels are active. Our aim here is to ascertain 
how their presence modifies the previous results obtained in the degenerate case, specifically 
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whether the system still lives in the phase where the U(l) symmetry is spontaneously broken. 
In this connection it helps, as we shall see, to derive an exact analytic expansion of the lowest 
eigenvalue of (pQ) in the inverse powers of the coupling constant g. This item has been lately 
pursued by Yuzbashyan et al. [13] in the framework of the Richardson equations. We also 
mention that this problem has been recently tackled in the path integral formalism as well [llj . 
We shall do the same here, but providing an expression for the expansion coefficients valid for 
any single particle level (s.p.l.) distribution and more transparently linked to the moments 
characterizing the latter. Moreover new results concerning g, namely the value of the coupling 
constant defining the domain where the expansion holds, will be presented. 

3.1. Comparison with the degenerate case 

We consider the ^-dependence of the ground state energy to assess how it is affected by the 
presence of the L single particle levels. For sake of illustration in this Subsection we confine 
ourselves to explore a model of equispaced single particle levels with unit energy distance and 
with unit pair degeneracy (Qj = 1). In Fig.[T]we display the Richardson solution for L = f2 = 12, 
for re=2, 6 (half filling), 8, 12 (full filling). The energy of the degenerate case (|12p is also shown 
as a reference. Moreover the g-behavior of the energy of each pair is displayed. 

It is clearly apparent in Fig. Q] the merging mechanism which allows the escaping of the pair 
energies from the grid of the single particle energies from where they originate at g = (see 
Refs. [7|2j ). In fact, considering the ground state, the pair energies start out, from the lowest 
bare single particle energies, being real at small g and then, two by two, merge at the energy 
of a single particle level (identified by the superscript i) for critical values of the coupling 5crit : 
clearly [n/2] of these exist. For g > g l cr[t the two pair energies associated with the index i become 
complex conjugate and, what is important, their common real part behaves almost linearly in g. 
Thus when g > g™it x an the pair energies become complex with real parts displaying a behavior 
close to linear in g and, remarkably, they add up to yield a downwards pointing straight line 
with the same slope as (fT2jh but upwardly shifted with respect to this one because of the 
contribution stemming from the single particle energies, which is obviously g- independent. One 
can accordingly argue that when L single particle levels are active, the system, at variance with 
the degenerate case, can live in different regimes. These in turn depend upon the filling of the 
levels. Here we consider the half filling situation where the system can live in three different 
regimes according to whether g < g™£ (normal fermionic regime, where the wave function of each 
pair has essentially only one component), g > g™it x (superfluid regime, where the wave function 
of each pair is spread out over all the L single particle unperturbed levels) and g^ < g < <?™it x 
(mixed regime). 

From the figure it also appears that the energy of the degenerate and non-degenerate case are 
markedly different only for g < <7™^ x . However this difference becomes less and less pronounced 
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as the ratio n/Q, approaches unity (full filling). Indeed for g < g™j^ x , but n sufficiently large, the 
non- linear behaviors of the exact E v (g) cancel out leading to a system's total energy linear in g 
as in the degenerate case. This occurrence relates to the absence of an excitation spectrum in 
the full filling situation in both the degenerate and non-degenerate cases. 

It is interesting to check whether the above findings agree with the Anderson criterion [12] for 
establishing how small a superconductor can be. This states that in the superconducting phase 
it must be 

1 > 1 ' (24) 
d being the average distance between the levels of the finite system and A the gap. Now by 
combining the above with (|23p it follows that a superconducting regime sets in for 



which indeed appears well fulfilled in our numerical analysis (which assumes d = 1, O = L and 
e given by (|27p ). Hence (|24p holds valid in our model. 

Finally a comment on the above recalled escaping mechanism, responsible for the setting up 
of a strongly collective ground state, is in order. 

Indeed the way the pair energies evade the grid changes its nature when the number of pairs 
substantially exceeds the half filling of the single particle levels. This is clearly apparent in the 
last panel of Fig. Q] (note that there a larger span of g is displayed), where the energies of the 
four upper pairs are seen to grow rather than to decrease with g. 

In fact the pair energies are differently affected by the pairing force. In particular the lowest 
lying pairs feel an attractive interaction whereas the upper lying ones feel a repulsive force. Since 
the sum of the pair energies is constrained to yield the degenerate result, it is not surprising 
that when n is large the upper lying pairs contribute to the total system's energy positively and 
not negatively. In other words the linear behavior in g of their energies has a positive (and not 
a negative) slope. 

3.2. The expansion of the ground state energy 

We now turn to derive an exact analytic expression for the 1/g expansion for the system 
ground state energy when L s.p.l. and n pairs are active. This we achieve along the lines 
illustrated in Ref . [13] , but we prefer to express the coefficients of the expansion in terms of the 
mean energy 




^Ve 2 + d 2 



(25) 




(26) 



and of the moments 




(27) 
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of the s.p.l. distribution (the first of these, namely when k = 2, 3 and 4, correspond to the 
variance, skewness and kurtosis). In the above O = X)u=i O^ i s the total pair degeneracy. This 
allows to better grasp the impact of s.p.l. on the spectrum of the system. 
We obtain 

00 fd\ j 

£ gs (n) = 2ne- gn(n-n + l) - d^aj (-) . (28) 

i=i 

In (I28p d is the average distance between the single particle levels and 

a i = fi E 2P+1 a S J™^ + 1 | P | 1 Je ?+1 - i m (,) J ■ (29) 



The coefficients a in (|29p are related through the following recurrence relations 

<+i = -^{| + (i-WEE <+U 4 + «(i - M E 2 fc a^ +1 E (J) , (30) 

which allows to evaluate the expansion (|28p at any given order starting from ag = nS^o- 

Eq. (|28() represents, when it converges, the Laurent expansion of £ gs around the simple pole 

it displays at g = oo. The first two terms of the expansion just correspond to the ground state 

energy of the degenerate case. The other terms thus account for the role of the s.p.l. distribution: 

indeed the n-th order in 1/g involves the moments of the s.p.l. distribution up to the order n + 1. 

Clearly the smaller g is, the more moments are needed to get a faithful representation of £ gs (n). 

This is supported by all the cases we have numerically explored. 

Keeping the first five terms only one gets (note that our term 1/g 2 differs from the one in 

Ref. \m) 



m^ 2 ) n(0 — n) ra(0 — n)(0 — 2n) 

£ gs {n) = 2ne - $n(0 -n + 1) - 2 x -- — -L + 2 x 



($0/2) 0-1 (5O/2) 2 (0-l)(0-2) 



2 n(0 - n) 



( 5 0/2) 3 (O- 1) 2 (0 -2)(0- 3) 
^O-SJ-ntO-n)^- 1 " 3 "- 4 )-" 



{ [0 2 (0 - 1) - n(0 - n)(50 - 6) 



(2) 2 



n-i - ■■ < 31) 

For one pair (|3ip yields back the result of Ref. [3] , whereas when all the levels are filled (n = O) 
it yields back (|12p . characterized by the linear ^-behavior displayed in Fig. [TJ Moreover, when 
the levels are equispaced only even moments, hence odd powers of 1/g, enter into the expansion. 

In Fig. [21 we compare the n-behavior of the ground state energy of the system as obtained 
in the Richardson's framework and from (|3ip . They are both displayed for a few values of g 
assuming L = 12 and the same set of equispaced single particle energies as before. From these 
results one can gauge the validity of the expansion (|28p even when only a few terms of the latter 
are kept. It is gratifying that few moments of the level distribution are enough to get a good 
representation of the true solution, providing g is not too close to the boundary of the domain 
of convergence of the expansion. 
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Accordingly we address the problem of establishing the domain of convergence of the expansion 
(|28p . To explore this issue we first consider the half filling situation in the following limiting 
cases: the macroscopic limit and when only one pair lives in two levels. In the former we assume 
that all the unperturbed levels have unit pair degeneracy and live in a finite energy band of 
width 2hujQ = 2dn = d£l, in such a way that when g — » the ratio g/d remains finite. We then 
obtain for the condensation energy 




3 



£•;;:>"' (n) 8, l An)--2nc-. -nh^ \± + ~ - ± (?) +-| (32) 



which coincides, when it converges, namely for 

5>! = l/», (33) 
with the Laurent <i/g-expansion of the BCS energy 

8%g$(n) = -nhLUocoth(-) (34) 



in accord with the long known result that the BCS energy cannot be computed perturbatively 
for small g/d. 

On the other hand the macroscopic limit of the Richardson solution has been explored in 
Ref. [7J, where it is shown that in this instance 

5 \ max 1 11345 J (9\ min 35 

rf/crit arcsinh(l) ' ' ' \dj crit 

Thus the domain of validity of the expansion (|32p does not directly relate to the critical values 
of the coupling constant: indeed although these correspond to branch points of the pair energies 
[6] these singularities cancel out in the expression for the system's global energy. 

Turning to the case of one pair living in two levels with unit degeneracy and having an energy 
distance d from Eq. (|3ip one gets 

e = - g -J 1 + l(i) 2 -l(i) 4 + ±(i)\..) (36) 



,9/ 8 \gj 16 \g 
which coincides, when it converges, namely for 

H = 1 > <37) 

with the well-known result 



£ = -g- • (38) 

Of course the range (|37p is not related to critical values of g which are now non existent. Note 
that this most simple case already shows that a g can exist without any g cri t. This occurrence 
goes in parallel with the findings of [2] where the same is seen to happen in a much more 
involved situation. 
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Thus for both systems above examined (n = 1 and n = oo at half filling) the function £ gs {g) 
displays a singularity (either a branch point or a pole) on the imaginary axis of the complex 
g-plane. By contrast the branch points associated with the g cr it all lie on the real axis. 

To gain insight on the actual value of | in a generic system we quote in Tables 1 and 2 the 
values of (4) mi f, 

V a i crit ' 



(d) 111 ^ an d °f d f° r a ^ ew vames °f n - The radius of convergence § of the 



expansion has been computed using the Cauchy-Hadamard criterion 



9 
d 



lim |a 7 '| 1//j 



(39) 



To assess the sensitivity to j of the above in Table 2 we report in the first column the results 
obtained with j = 21 and in the second column those obtained (at least in a few cases) with 
j = 199. 

From our results it follows that the radius of convergence •? of the expansion (|28|) for any n 
appears to occur, for the half-filling case and for equispaced s.p.l. of unit pair degeneracy, in 
the domain ranging from (f)™™ to (f )™^ x - This interval increases with the pair number n, as 
it can be inferred from Table 1, where (§)™^ x ((f ) r -f ) is seen to increase (decrease) with n. 

We conclude that | 

a) varies in the range - < % < 1, for any n, 

b) decreases with n, 

c) lies, for a given n, in the range (f)™™ < § < (f ) 



max 
crit 



Note that ^ , even for n = 500, is still far from 1/ir : hence the macroscopic limit of the pairing 
Hamiltonian appears to be reached very slowly. 

We now abandon the assumption of s.p.l. of unit pair degeneracy and address the question of 
the impact on | of the pair degeneracy of the s.p.l.. In the simple case of one pair in two levels 
the generalization of (|38p is easily found to read 

fii + ^2 



£ = (ei + e 2 ) - g- 



which yields 



ill + ^2 



+ d 2 + gd(Q, 2 - Hi 



n 1 -n 2 ± 2i^/^k^h 



(iii + n 2 ) 2 



(40) 



(41) 



being d = ei — e 2 , from where it is seen that the singularity moves away from the imaginary 
axis when Hi ^ il 2 . 

Sticking to the case Hi = il 2 = il one gets from (|41[) 
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n 


t g \ min 
Vd/crit 


/ g \ max 
Vd/crit 


2 


0.66 




4 


0.46 


0.82 


6 


0.40 


0.90 


8 


0.35 


0.94 


10 


0.32 


0.96 


20 


0.27 


1.03 


OO 





1.13 



Table 1 

For a few values of the pair number n the largest and smallest critical values of the coupling 
constant are quoted at half filling. The number of critical points is n/2 and clearly for n = 2 
only one of them exists. 



n 


|«2l| 1/21 


U 1 1/199 
|Ol99| ' 


1 


0.82 


0.96 


2 


0.69 




50 


0.41 




100 


0.39 


0.65 


500 


0.36 


0.64 


oo 


0.31 


0.318 ~ l/vr 



Table 2 

The radius of convergence of the expansion (|28p for various pair numbers n at half-filling and 
equispaced single particle levels of unit pair degeneracy. 

and from (|4"Dj) (choosing e\ + ei = 0) 

£ = - g n - ^d 2 + g 2 n 2 (43) 

which coincides with (|38p when = 1. Thus the larger 0, is, the closer to the real axis the 
branch point is. This occurrence has a bearance on the physics occurring on the real axis: in 
fact not only the energy but also the pair wave function is strongly modified by O, in the sense 
that the larger O is the sooner in g the two components of ([5]) become equally weighted. 

In the limit — > oo where the pair becomes a true boson the singularity lies on the real axis. 

To assess the impact of the s.p.l. filling on g/d we quote in Table 3 the numerical values of 
the latter in the case of two levels with equal pair degeneracy f2i = O2 = 100 for a few values of 
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n 


U 11/199 
|«199| 


1 


0.0096 


50 


0.0077 


100 


0.0068 


150 


0.0077 


199 


0.0096 


200 






Table 3 

Values of | in the case of two levels with equal pair degeneracy fii = i\ = 100 for a few values 
of n. 



n. From Table 3 one sees that g/d reaches its minimum at half filling being symmetric around 
the minimum. 

Finally we give a striking example illustrating the crucial role played by the s.p.l. degeneracy 
on ~g/d. For this purpose we consider the two following cases 

a) L = 2, Qi = Q 2 = 100, Q Tot = 200, d = 1, n = 100 
and 

b) L = 200, fii = 2 = ■ • • = 1, ^Tot = 200, d = 1/199, n = 100. 

Then from the explicit computation it turns out that in a) g/d = 0.0068 (see Table 3) whereas 
in b) ~g/d = 0.65. Thus even in situations where analytical results are hard to achieve, the 
numerical analysis suggests that g/d oc being O, the degeneracy of each s.p.l.. 

In concluding this Section we return to the problem of the thermodynamic limit previously 
mentioned pointing out that two quite different variational solutions, namely the Richardson (JH) 
and the BCS ([3]) projected on a given particle number M which reads 

yield the same ground state energy. Whether this remains true for other observables as well is 
an issue worth to be further explored. 

4. The spectrum of the pairing Hamiltonian 

In this Section we shortly examine the spectrum of the pairing Hamiltonian in the non- 
degenerate situation where ([!]) predicts new excited states (referred to as Gaudin excitations) 
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beyond those associated with the breaking of pairs (referred to as seniority excitations), which 
are the only ones present in the degenerate case. 

These new kind of excited states correspond to the raising of pairs into higher-lying s.p.l. and 
are classified in terms of the Gaudin's number Nq [15], which yields the number of pair energies 
staying finite when g goes to infinity because they remain trapped in the grid of the s.p.l.. For 
fixed Nq (of course < Nq < n), the number of such excited states is (^ ) — (jy ])■ 

As an illustration of these excitations we consider the simple example of 3 pairs in 5 equispaced 
levels of unit pair degeneracy. In this case there are 4 states with Nq = 1 and 5 with Nq = 2 
(note that Nq = 3 is not allowed). In Fig. [3] we display the pair energies versus g for the ground 
state (panel a) and for the Nq = 1 states (panels b-e) together with the associated total energies 
(panel /). Note that while for the ground state and for each of the first three excited states only 
one (7 cr it exists, for the fourth Nq = 1 excited state (panel e) a complicated escaping mechanism 
takes place, giving rise to three g cr it- In the last panel, where the total energy of the above 
states is plotted, it clearly appears that, for g high enough, the states tend to group according 
to their Gaudin number in accord with the finding of |13| . 

Turning to the seniority excitations of ([1]), their number, for levels with unit pair degeneracies, 
is („ ) ( ~_ 2 T) (of course it must be < s < min{n, $7 — n}). 

Again for illustration we consider the very simple set of configurations associated with 2 pairs 
living in 4 equispaced single particle levels with unit energy distance and pair degeneracy. In 
this case 13 seniority excited states occur, 12 with s = 1 and 1 with s = 2. Of the former, 
owing to the degeneracy entailed by the symmetry of ([1]) , only 6 show up in Fig. [J^i, where the 
g-behavior of the seniority eigenvalues is displayed. These can be analytically computed since 
they are obtained as a solution of a two-levels problem, the other two levels being blocked by the 
broken pair. Assuming, e.g., that the blocked single particle levels correspond to the energies e\ 
and e2 (the other cases obtain with trivial permutations) they read 

£{n = 2, s = 1) = ei + e 2 + e 3 + e 4 - g ± yfg~ 2 + (e 3 - e 4 ) 2 , (45) 

with obvious meaning of the symbols. The ± sign accounts for the splitting of the s = 1 states 
into two families. It turns out that the energy of the lower lying family comes pretty close to 
the energy of the Nq = 1, s = state (compare Fig. Sb), whereas the energy of the higher lying 
family comes close to the one of the Nq = 2 state. Moreover these are almost degenerate with 
the s = 2 state, corresponding to a configuration where all the four levels are blocked: hence 
their energy is simply given by the sum of the unperturbed single particle energies. 

Generalizing the above example we argue that the excited states of the pairing Hamiltonian 
are identified by two quantum numbers, namely Nq and s. However the associated eigenvalues 
group, as g becomes large, into families characterized by the number iV^ of pairs blind to the 
pairing interaction: whether they are so because they are broken or because they are trapped it 
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does not matter for g sufficiently large. 

For the energies of these states, inspired by ( (|13p ). we accordingly surmise the following 
approximate expression 

£(n,N B ) = S 9 . s .(n) ~ W^b) + 2N se , (46) 



which we test against the exact solution in the N B = 1 case. Confining ourselves to (|3ip the 
above becomes 

2" 



£ (n, N B = 1) - £(n, N B = 0) = 2 



g n m(4) 



(2) 



2 (pfi/2) ( 5 n/2)2 ( 5 fi/2)S 



(47) 



whose first term, significantly, is equal to 2 A. 

Notwithstanding the shortcomings of formula (|46l) yet we have found, as seen in Fig. 0b, that 
it provides a good account for the exact results, except at low g namely in the domain of the 
critical values of the coupling (in our example these occur at 0.66, 1, 2 and 3 in units of the 
spacing between the single particle levels). In addition (|47p shows that to leading order, which 
is the accuracy of formula (|47p . the Gaudin and the seniority excitations are both characterized 
by the same gap £lg. 

5. Conclusions 

The problem of the pairing Hamiltonian, extensively treated in a number of papers, is still the 
object of many investigations. Such an Hamiltonian employs a drastically simplified interaction 
and, as a consequence, cannot be viewed as a realistic one. It yields, for example, a gap too 
large by a factor (4e) 1 ^ 3 with respect to the accurate predictions of Gorkov [16] and Heiselberg 
|17j . Yet the fact of being exactly solvable makes (pQ) interesting and helpful for the researches 
in a variety of fields. To mention a couple of these, we remind the microscopic derivation of the 
Arima-Iachello model and the physics of small superconducting metallic grains, which has lately 
attracted much interest. 

In this work we have addressed a few aspects of the problem which, in our view, still deserve 
further studies. First we have proved in the path integral framework and in the degenerate case 
the superfluidity of the pairing Hamiltonian solution in a "finite" system. This item was already 
lately treated by us pQ. Here, however, we have added the important finding that the saddle 
point of the effective action derived in [1] is exactly fixed by the BCS gap. 

Then we have compared the (/-behavior of the ground state energy in the degenerate case 
with the situation when many single particle levels are active. Although this comparison has 
been limited to levels with unit pair degeneracy it has allowed us to grasp the impact of the 
s.p.l. distribution on the system's energy and to get, at least in the half filling situation, an 
orientation on the regimes in which the system lives. 
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Next we have derived an expansion of the system ground state energy £ g . s {g) in the inverse of 
the coupling constant (strong coupling expansion) . The formula we have obtained generalizes the 
one of Ref. [13] since it is expressed in terms of the statistical moments of the s.p.l. distribution 
and is thus valid for any distribution. This appears to be useful for some recent researches on 
the density functional approach to the pairing Hamiltonian [T8] , but it is also helpful for a deeper 
understanding of the dynamics of a system ruled by (JTJ) . Indeed the investigation of the domain 
of validity of the expansion has lead us to conclude that such a domain is set by a singularity 
of the function £{g) located somewhere in the complex g-plane. Although we have not been 
able to find a general expression for this singularity (indeed it depends from the specificity of 
each case) we have numerically found where its modulus g lies at least for s.p.l. with unit pair 
degeneracy and how its location is affected by the distribution of the levels (in particular by 
their degeneracy) defining the space on which the Hamiltonian (fTJ) acts. This is of relevance 
because we conjecture that the physics is more or less affected by the singularity depending 
upon its proximity to the positive real axis, as it has been hinted by the simple n = 1, L = 2 
case: accordingly we surmise that if the singularity is located remotely from the real axis one 
deals with a system which is predominantly fermionic, if it is closer to the real axis (as it may 
be the case when the s.p.l. have large pair degeneracies) one deals with a system in which a 
process of bosonization is taking place. 

Finally we have discussed the excitation spectrum of the pairing Hamiltonian. We have found 
that the Gaudin and seniority excitations are close to each other, at least for small s and Ng. We 
have proposed an approximate analytic expression for these excitations that appears to account 
satisfactorily for the exact result. Importantly, both the Gaudin and the seniority excitations 
display the same gap in the leading order expressed by our approximation. 
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Figure 1. The real part of the exact ground state solutions E v (thick solid curves) of the 
Richardson's equations (6) for a system of n=2,6,8 and 12 pairs living in a set of L=12 equispaced 
non-degenerate (Qk = 1) single particle levels with energies e^ = k (k = 0, • ■ ■ 11) are displayed 
versus the strength g together with their sum £ (n) (thin solid curves); the dotted curves represent 
the shifted total energy £(n) — 2J2k=i e k and the dashed curves the energy of n pairs living on one 
level of energy e = and degeneracy L. The horizontal dashed lines represent the unperturbed 
solutions and the cuspids are due to numerical imprecision close to the critical points. These 
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Figure 2. Left panel: the exact (solid curves) ground state energy for the same model as in 
Fig. [1] is displayed versus the pair number n for g=l (upper curve), 2, 4 and 6 (lower curve), 
respectively; the dashed curves (visible only in the g=l case) correspond to the approximate 
solution of Eq. (29). Right panel: the same curves divided by g; for very large g the curves tend 
to the parabola —n(L — n + 1). 
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Figure 3. The real part of the exact solutions E u of the Richardson's equations (6) for a system 
of n=3 pairs living in a set of L=5 equispaced non-degenerate (Qk = 1) single particle levels 
with energies ej. = k {k = 0, • • • 4) are displayed versus the strength g. Panel a: ground state; 
Panels b-e: excited states with Gaudin's number Nq = 1; panel /: total energy of the ground 
state (lower curve) and of the Nq = 1 states (upper curves). 
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Figure 4. a) Exact energies of the ground state (solid), seniority s = 1 (dashed) and s = 2 
(dotted) states for a system of 2 pairs in 4 equispaced single particle levels of energies 0, 1, 
2 and 3, respectively, and pair degeneracy 1. The s = 2 state trivially corresponds to E = 
ei + e2 + e3 + e4 = 6. b) Exact energies of the s = states with Nq = (solid), Nq = 1 
(dashed) and Nq = 2 (dotted), c) Comparison between the Nq = 1, s = 1 states (dashed) and 
the Nq = 1, s = (solid) excited states. The crosses correspond to the approximate solution 
(37) for N r = l. 



